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SUMMARY 


In previous work by Pratt and Whitney and NASA-Lewis, a rotating microphone data acquisition/ 
reduction system was developed to resolve acoustic pressures in an engine inlet plane into constitu- 
ent harmonics of blade passing frequency and circumferential and radial mode orders. Although 
this system met its design objectives, it had no means of distinguishing forward (upstream-going) 
waves from waves that might be reflected downstream from the engine inlet. 

This report describes a procedure for enhancing the use of the basic rotating microphone system so 
as to determine the forward propagating mode components of the acoustic field in the inlet duct at 
the microphone plane for more accurate prediction of far-field radiation patterns. In addition, a 
modification was developed to obtain, from the same microphone readings, the forward acoustic 
modes generated at the fan face, which is generally some distance downstream of the microphone 
plane. Both these procedures employ computer-simulated calibrations of sound propagation in the 
inlet duct, based upon a finite element radiation code. 

These enhancement procedures were applied to previously obtained rotating microphone data for 
the 17-inch ADP fan. The forward mode components at the microphone plane were obtained and 
were used to compute corresponding far-field directivities. In-duct mode amplitudes and far field 
directivities were compared for the basic and enhanced systems. For the four cases evaluated, the 
enhanced data reduction scheme made only a small difference in predicted noise levels. The reason 
for this negative result is that, while some modes were near cutoff, the microphone array was not 
located near a pressure extreme in the standing wave pattern. 

To demonstrate that the enhanced data reduction scheme is needed under conditions where a mode 
might be near cutoff, an analytical simulation was performed using an acoustic field inside and out- 
side an inlet generated by a finite element acoustic code. The basic and enhanced methods for proj- 
ecting internal measurements to the far field were applied, assuming various locations for the rotat- 
ing microphone plane. Depending on mic plane location, far field predictions using the basic data 
reduction scheme varied over a range of 12 dB. When the enhanced scheme was applied, the far 
field directivity patterns collapsed into a ± 0.5 dB band. The success of this demonstration indicated 
that the reflected wave effect can introduce large error into estimates of the outgoing modes from 
the fan. It also indicated that the new scheme can be very effective in eliminating the reflected wave 
contamination. 

The second main task of the program involved finding the forward wave modes generated at the fan 
face in terms of the same total radial mode structure measured at the microphone plane (nearer to 
the inlet). This was not as simple a matter as before; initial calculations gave obviously wrong re- 
sults. To resolve this problem, a set of computer-simulated sensitivity studies was made that showed 
the effects of small measurement errors. These indicated that, for the spacing between fan and mi- 
crophone planes used in the ADP tests, the enhancement procedure could be applied successfully 
only to the propagating radial modes. With this limitation, satisfactory predictions of the fan modes 
were obtained using the test data. 

The sensitivity studies also led to recommendations for future design and use of the rotating micro- 
phone system on other fan configurations. 



INTRODUCTION 


The rotating microphone system (Reference 1) was designed to provide detailed information on the 
circumferential-radial mode structure of the acoustic field in inlets of turbofan engines and to obtain 
this information using a number of microphone elements that is dramatically smaller than would be 
needed with fixed microphone arrays. 

The concept was applied to a rotating rake which provided 9 radial microphone locations. The mi- 
crophone signals were frequency— analyzed with a very high resolution Fourier transform. Each fre- 
quency thus obtained could be uniquely associated with the spinning mode index (m); this separation 
is made possible by the (low) rotation speed of the microphone array (Reference 1). For each spin- 
ning mode sensed by the microphones, the associated radial mode structure is obtained by least 
squares fitting the 9 complex microphone pressures to the first 4 associated radial mode (Bessel) 
functions. 

In 1991 rotating microphone tests were run on the 17-inch diameter Advanced Ducted Propulsor 
(ADP) rig in a wind tunnel at Lewis Research Center. Descriptions of the test procedures and rotat- 
ing microphone usage are given in References 2-4. 

These tests provided the circumferential and radial mode structure of the total acoustic pressure in 
the inlet at the plane of microphone location, and were also used as a basis for computing the 
associated far field radiation directivity patttems for selected modes. However, the pressure mea- 
sured by the rotating rake consists of two components: The primary component is the forward going 
wave generated by the fan, p + , but there is also present a backward wave reflected from the mouth 
of the inlet, p~. Only the p + forward wave generates the outside radiation field, and use of the total 
measured pressure, p = p + + p", will introduce error if p - is significant. 

For the measured cases which were studied, far field differences were not large. In an effort to deter- 
mine if there are conditions for which large differences can exist between the basic and enhanced 
systems, an analytical simulation was run with the results shown in Figure 1. The geometry used 
for the simulation was the ADP long inlet. At a speed of 7500 rpm. Blade Passing Frequency propa- 
gates in the 6,0 mode which is just slightly cut on. (Cutoff ratio £ = 1.12 at the fan face and 1.07 
at plane M.) All other modes are well cut off and can be neglected. The noise field within the duct 
was predicted assuming a unit amplitude input at the fan face for the single propagating mode and 
is shown in the top panel of Figure 1. Following this, the noise field was interpolated to obtain pres- 
sure data as though measured by a microphone array for the three other axial locations indicated in 
the figure. The enhanced system was then used to determine the p + levels at each array location. 
Far field estimates based on both p and p + at the source and the three other axial locations were com- 
puted and are shown in Figure 1 . Differences of about 12 dB can be seen for the data estimated from 
p (center panel) while the predictions based on p + (bottom panel) agree to within about 0.5 dB . Near 
a pressure maximum such as location H, the p levels can approach p + + 6 dB. Near a node such as 
location G, the total p can become very small compared to p + , resulting in a very large change in 
predicted far field levels. 

The remainder of this report presents a derivation of the equations for the enhanced system plus ap- 
plications to rotating mic data from the 17 inch ADP model. 
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TASK DESCRIPTION - ENHANCED ROTATING MICROPHONE METHODS 


The procedure for enhancing measured microphone outputs involves using Eversman’s radiation 
code (Reference 5). By using finite element numerical methods, this computes the acoustic field in 
the duct and the far field that results from known forward wave inputs. Inside the inlet, this field 
is considered to be the sum of both forward-going and reflected waves. 

The task was divided into two parts for convenience: 

Application 1 - Improved method to predict far field noise from rotating microphone data. 

This application has been described in the Introduction. It consists of devising a way for extracting 
the forward wave component, p + , from the total measured pressure at the rotating microphone plane 
and using this quantity to obtain an improved prediction of the far field. 

Application 2 - Source diagnosis method 

It turns out that the basic methodology involved in Application 1 can be modified in a fairly obvious 
way to attack another important problem in fan noise measurement and interpretation. This is the 
matter of determining what forward modes the fan is generating as a source. This information is 
needed to verify theories about fan noise generating processes and their prediction. If the rotating 
microphone could be placed close to the fan leading edge plane, the methods of Application 1 would 
suffice. However, mechanical constraints prevent this extreme positioning. Instead, a way must be 
found to use microphone data measured in a plane nearer to the inlet to infer the source characteris- 
tics. 
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ANALYSIS 


Basic Relations 

In any plane normal to the fan duct axis the periodic acoustic pressure is represented here for simplic- 
ity as the real part of 


00 00 

P(0, r, t) = X X expi(m0 - cot) ( 1 ) 

m= - oo [i = 0 

where ^m^r) = J m (k m(1 r) + Q m(l Y m (k mti r) , and the eigenvalues k mfl and Q m(1 are evaluated for 
the appropriate hub and tip boundary conditions at that plane. 

Operation of the rotating microphone system automatically separates frequency (w = BPF, 2BPF, 
etc.) and circumferential mode wavenumber, m. Then for any m and ( 0 , Equation ( 1 ) reduces to 


P(r) = £ C m ^ m| x(r) (2) 

fi=0 

Values of the complex coefficients, C m(l , are obtained from measured values of P(r) taken from the 
microphone outputs at a number of radial locations (9 in this program). A least squares fit of Equa- 
tion ( 2 ) to the 9 values of P(r) serves to give values of C m(i . While, in principle, it is possible to 
obtain the first 9 coefficients from 9 points, it is well known , e.g. (Reference 6), that slight measure- 
ment and/or computation errors can render the results useless. By fitting Equation ( 2 ) to only about 
one-half the maximum number of modes, more reliable results are obtained. For this reason a fit 
to the first 4 radial modes has been used in this program. If, in other applications, a larger number 
of modes are known to be propagating, provision must be made in advance to employ a suitably larg- 
er number of microphone locations. 

In the following work it will simplify and clarify matters if we make the following change in nota- 
tion: With the understanding that m and oa are known and fixed, a radial mode can be identified by 
the coefficient C m(i : this can further be abbreviated by replacing C m(i by the symbol p^. 

Complex coefficient of fx-th radial mode = p^ ( = C m(i ) ( 3 ) 

If several radial modes are present in a transverse plane, the collection can be represented simply 
by a vector or column matrix: 


For example, with 3 modes 


Pi 

P2 

P3 
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P 


(4) 



p represents the measured total pressure field in a transverse plane and is extracted from the data as 
has previously been described. It is the result of waves traveling upstream from the fan and reflected 
waves traveling downstream from the inlet, or p + and p“ waves: 


that is 


for 3 modes. 


p = p + + p , 


Pi" 

= 

Pl + 

+ 

Pi- 

P2 


P2 + 


Pi" 

P3 


. P 3 + _ 


P3 _ 


(5) 


The general objective of this investigation is to separate the measured total mode pressure, p, into 
its forward and backward wave components, p + and p _ . 

The key to solution lies in “calibrating” the duct with a sequence of unit forward wave inputs, p + , 
and finding the resulting outputs, p. These calibrations will be made by use of the radiation code 
as explained below. 

The result of the calibration may be expressed as the linear relation 

p = Dp + , (6) 

where elements of the square matrix, D, are now considered known. 

The solution is obtained by inverting ( 6 ) to give the desired forward waves, p + , in terms of the mea- 
sured waves, p: 

p+ = D _1 p. ( 7 ) 

The reflected waves can be obtained from 

p- = p-p+ = (I - D -1 ) p. (8) 

Eversman’s code is applicable since it gives p everywhere in the acoustic field for specified values 
of p + input at the boundary of his computational domain, a plane in the inlet. 

Details of the calibration procedure and its use in processsing measured data, p, will now be pres- 
ented, first for Application 1. 
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Analysis for Application 1 — Improved Method to Predict Far Field Noise from Rotating 

Microphone Data : 

In previous work the measured microphone mode structure has been used as input to the radiation 
code to compute far field directivity. This application requires that the total mode structure at the 
microphone be resolved into forward and backward waves, and that only the forward component 
be used to compute far field. 


To signify that the microphone plane specifically is involved, subscripts M are appended to Equa- 
tions (5), (6), and (7), giving: 


Pm = Pm + Pm 

(9) 

Pm = d mm Pm 

(10) 

Pm = d mm Pm 

(ID 


These equations are illustrated in Figure 2. 


To find the elements of Dmm> the radiation code is set up with the origin of the mesh system located 
at the microphone plane. This mesh system is shown in Figure 3. The code computes values of the 
(total) pressure, p, in the inlet for any specified p^J mode structure. 


To “calibrate” the duct, a sequence of unit p^, modes is applied at the mesh origin and the resulting 
total pressures are computed at points in the microphone plane. These pressures are converted into 
radial mode constituents by least squares methods. (43 mesh points were used in the microphone 
plane, so that the least squares procedure used in computer simulated calibration will determine the 
first few modes to a high degree of accuracy.) 

Each unit mode vector, p^, will produce output mode values, p M , which constitute the elements in 
the corresponding columns of the D-matrix. Figure 4 illustrates this calibration procedure, using 
a 3-mode system for convenience. 

With the D-matrix in hand, it is then possible to extract the forward wave components from rotating 
microphone data taken in tests on the 17-inch ADP rig. Details of these procedures and results are 
described in a later section of this report. 


Analysis for Application 2 - Source Diagnosis Method 

In this application the same microphone readings, p M , are to be used to infer the forward mode struc- 
ture generated at the fan face . This information will be of use in validating theoretical models for 
noise generation by rotor-stator interaction. 
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The relation between microphone measured modes, p M , and the fan source forward modes, denoted 
here by pp , will have a linear structure given by the (different) matrix equation 


Pm ~ d mf Pf 


( 12 ) 


Figure 5 shows this relation for Application 2. 

A new calibration sequence of unit pp modes must be used to establish values of the matrix 
elements by setting up the radiation code mesh origin back at the fan face. This mesh is shown in 
Figure 6. 

The calibration procedure for Application 2 is shown in Figure 7. (This procedure may be easier 
to visualize, since the forward wave input, pp , and the total pressure output, p M , are located at 2 
distinct places in the duct.) 

With the elements of Dj^p in hand, the desired fan mode structure is obtained by inversion: 



( 13 ) 


Procedures and results involved in Application 2 are described in the next section. 



PROCEDURES AND RESULTS 


The enhancement procedures described in the Analysis section were applied to the geometry of the 
17-inch ADP rig incorporating the so-called “medium inlet” or “mid-inlet”. See Figure 8. Data 
from the rotating microphone array were extracted at fan speeds of 9,600 rpm, 11,400 rpm, and 
12,000 rpm. For the configuration incorporating a 22-vane stator with the 16-blade rotor, the inter- 
action tone has a circumferential mode m=-6 (=16 - 22) for the fundamental blade-passage order, 
BPF. Additionally, at 9,600 rpm 2BPF was analyzed for m=10 (=32-22). Radial mode structures 
were obtained by NASA LeRC, using a least-squares fit of 4 radial modes to the data at 9 radial 
microphone locations. In Figures 9 through 12 the demarcations between propagating and decaying 
modes are indicated. 


Application 1 

This application consists of extracting the forward mode component vector, p^, from the measured 
(total) microphone mode structure, (p M ), f° r the test conditions described above. This pjJ structure 
was then used to compute the far field and to compare results with those obtained using p M as a 
source. 

The enhancement procedure involves obtaining the matrix in Equation 10 by computer simu- 
lated calibrations and applying it to the data, p M . The calibrations were made by using the radiation 
code with the source (input boundary) located in the microphone plane. The finite element grid is 
shown in Figure 3. 

Figures 9, 10, and 11 show the amplitudes of the 4 components of pjJ, together with p M , and the 
far field directivity for the sum of the complex components. These figures are for BPF at 9,600, 
1 1,400, and 12,000 fan rpm. Figure 12 presents corresponding information for 2BPF at 9,600 rpm. 

Some observations may be made about these 4 figures. For the conditions presented, the differences 
between mode amplitudes, p^ and p M , are less than 4 dB. In the majority of cases the magnitude 
of p M is greater than the forward wave, p^. 

In Figure 9 the directivity curves are indistinguishable. For this case there is only one propagating 
mode which is well cut on at this speed so that the reflected wave is practically zero. Since the for- 
ward mode component, p^, is nearly equal to the total mode amplitude, p M , this result is not surpris- 
ing. 

More significant differences can be seen in Figure 10, where two radial orders propagate and the 
p M and pj£ mode amplitudes for the /t=l radial order differ by about 3 dB at radiation angles be- 
tween 60 degrees and 90 degrees. At this speed the /z=l mode has a cutoff ratio of about 1.07 so 
that some reflection of this mode might be expected. At about 40 degrees where the major lobe peaks 
for the n = 0 mode, the p^ directivity pattern is slightly lower than the pattern generated using p M . 
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The difference is similar to the difference between the p M and the mode amplitudes for the /U=0 

mode. At larger angles, where higher order radial modes have their maximum amplitudes, the differ- 
ence between the directivity curves is similar to the difference in the fi=l mode amplitudes. 


At 12,000 rpm. Figure 11, the propagating modes are similar in amplitude so that little difference 
in directivities would be expected and this is seen to be the case. A larger difference is seen betweeen 
the n=2 mode amplitudes, but, since these are cut off, no differences would be expected to occur 
in the far field. 


Three radial modes propagate for the case shown in Figure 12. The ju=2 mode has a cutoff ratio of 
about 1.06 and appears to have a significant reflected component. This is likely the cause of the 
directivity curve differences above 50 degrees. Below 50 degrees, three modes with similar ampli- 
tudes combine and their relative phasing is probably responsible for the larger directivity levels for 
the p^ modes even though all three propagating mode magnitudes are lower. 


From these sample results it seems that far field directivity can be significantly changed if one of 
the modes is near cutoff. If few modes are present, differences will probably be observed at angles 
larger than that of the main lobe of the lowest order mode. Also, for airplane noise calculations, 
constant radius data is extrapolated to constant sideline distances, giving the results at larger angles 
somewhat more importance than is indicated here. As the number of modes increases, however, the 
effect of substituting p^ mode amplitudes for p M mode amplitudes will probably be harder to ob- 
serve since the directivity curves will be dominated by propagating modes and the modes near cutoff 
will not provide a large contribution to the total sound field. 


The pjJ mode amplitudes are, in principle, more appropriate for use in computing far field directivi- 
ties because these are the modes used in the formulation of the radiation code. For cases containing 
well cut on modes, the use of p M mode amplitudes will do an adequate job of estimating far field 
directivities with a significant saving of the computational effort involved in determining the pjJ 
mode amplitudes. For cases with relatively few modes propagating with one near cutoff, determina- 
tion of the appropriate p^ mode amplitude should be used. 


One assumption has been made in the calibration procedure for obtaining the matrix elements: 

As presently written, the radiation code takes the coordinate system origin (x/r = 0) as common to 
both the inlet duct flow and the acoustic field. Thus, in calibrating with the sound source in the mi- 
crophone plane, which is necessary, the axial flow in this plane is automatically constrained by the 
code procedure to be constant from hub to tip. This uniform distribution is an artifact of the code. 

Overcoming this limitation would require revising the code so that the flow origin is taken at the 
fan face, and then imposing the calibrating acoustic modes at the microphone. Feasibility of making 
such a change in the radiation code might be investigated in the future, but whether the results would 
be sufficiently different from the present ones can not be predicted in advance. 
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Application 2 


As described in the Analysis section, the objective of this application is to determine the modal out- 
put at the fan face, pp" , using the same microphone data, p M , that has been presented. It should be 
mentioned that the numerical procedure followed to obtain the elements of D^p differs in one signif- 
icant aspect from that used in Application 1: 

In Application 1 the source of the p + waves coincides with the radial plane containing the micro- 
phone pressures, p^. This plane, shown in Figure 3, is the left-hand boundary of the finite element 
mesh, and there were 43 mesh points extending radially from the inner radius of the duct to outer 
radius, at which the radiation code computes the total pressure, p M . 

For Application 2, however, the mesh starts at the fan face, downstream from the microphone plane. 
Figure 6 shows the left boundary of the mesh being in a radial plane. But it is also seen that the micro- 
phone plane cuts through the mesh cells in a complicated way, due to the mesh curvature. Since 
pressures are required at points in the radial microphone plane, and the radiation code gives the pres- 
sures along cell boundaries, a 2-dimensional interpolation is needed. A substantial effort was in- 
volved in devising and applying a suitable procedure. For this procedure a second degree polyno- 
mial in variables x and r with complex coefficients was chosen to represent the pressure in the 
neighborhood of each point in the radial microphone plane for which the pressure must be interpo- 
lated: 


p(x,r) = ai + a 2 x + a 3 x 2 + a 4 r + asx r + a^ r 2 

Values of the 6 complex coefficients were obtained by fitting this function to the 6 points on the mesh 
closest to the interpolation point, (x,r). This required solving 6 simultaneous equations for each of 
the interpolation points. 

Note: In the above descriptions, statements are made that pressures computed by the radiation code 
lie along the mesh lines. This was done to simplify the presentation for a casual reader. However, 
anyone intending to use the enhanced microphone system, should be aware that the pressures are 
really jnql calculated along mesh lines, but rather, at so-called “Gauss points”, which happen to lie 
inside of each cell. The subject of Gauss points is handled in Reference 5) 

Calibrations were then made to obtain elements for the same ADP conditions used in Applica- 

tion 1: Blade passage frequency at 9,600 ipm, 11,400 rpm, 12,000 rpm, and 2BPF at 9,600 rpm. 


The fan-generated modes, pp , were computed by Equation 13 from the measured microphone val- 
ues, p M . Table 1 lists p M and pp" for the first 4 radial modes of BPF at the 3 rig speeds. 

The listed numbers are the decibel values, re. one Pascal, of the mode amplitudes, and the enormous 
differences between p M and pp" for the propagating modes are an indication of a serious problem 
with the analysis system. In Application 1, p M and its component, p^, were found to be quite simi- 
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lar, so for a propagating mode, the fan pressure pp should not be nearly as high as listed in Table 
1. To verify this contention, calculations were made, using the radiation code, of the duct pressure, 
p, at a sequence of axial positions, x, upstream of the fan resulting from unit values of fan-generated 
modes, pp . 

Samples of these calculated results are plotted in Figures 13 and 14 for 9,600 rpm and 12,000 rpm, 
respectively. In both Figures 13 and 14 the designations a, b, c, and d indicate that the single unit 
fan forward mode applied to the system was the first, second, third, and fourth, successively. Incre- 
ments of 0.05 units of x (normalized by duct radius) were taken between x/r = 0 and 0.8. (The micro- 
phone plane in the ADP rig tests was at x/r = 0.6.) 

These figures clearly show that the levels of propagating modes vary only moderately from the ori- 
gin at x/r = 0 to microphone location (x/r = 0.6). (Modes below cutoff are shown to decay exponen- 
tially, with higher modes decaying more rapidly.) 

Consequently, the values of pp in Table 1 inferred from measured microphone data, p M , are ob- 
viously wrong. 

To clarify matters a sensitivity study was conducted that indicated the effect upon pp of microphone 
errors in the measured test data, p M . Computer simulation using known inputs, pp , to the finite 
element code was used to generate presumably accurate values of the total microphone pressure, p M . 
These outputs were then perturbed by amounts ranging from 0.1 dB to 0.5 dB to simulate errors in 
p M due to data acquisition and processing in an actual test. These perturbed values of p M were then 
used in Equation 13 to obtain the source fan modes, pp . Differences between these pp and the ini- 
tially known pp values provided indications of the sensitivity of the process to microphone mea- 
surement error. 

Details of these sensitivity studies arc presented in Appendix B. While, in the rig tests, the micro- 
phone plane was located at a distance x/r of 0.6, the computer simulated runs were made for a series 
of hypothetical microphone locations of x/r = 0 to 0.8 in steps of 0.05 to reveal detailed behavior 
of the system as the microphone plane successively moved farther along the axis from the fan. 

The fan-to-microphone separation turned out to be a critical parameter. For separations of less than 
x/r — 0 . 2 , good results were obtained using the first 4 radial modes, and these improved as x/r ap- 
proached zero. (This characteristic showed that the procedures involved in Application 1 were reli- 
able since the source-to-microphone distance used in computer calibrations was essentially zero 
there.) 

For distances greater than x/r = 0.2 results progressively degenerated; the computed values of pp 
differed more and more from the known input values. For the actual test microphone location at 
x/r = 0.6, the computed values of Pp had no recognizable relation to the initial known fan mode 
inputs. 
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On the basis of a conjecture that low levels of decaying modes at large microphone distances were 
responsible for the problem, the procedure was repeated with an important modification: Instead 
of using 4 radial modes to fit microphone point data it was decided to resolve the point data into only 
the propagating modes ( or, into those plus a weakly decaying mode, if present.) 

This modification worked. The perturbed microphone data gave back the known inputs satisfactori- 
ly for x/r = 0.6 and beyond, as seen in Appendix B. 

On the basis of the sensitivity studies, it is possible to select a least squares fit to give reasonable 
results for fan noise prediction, pp" , from the actual microphone test data, where x/r = 0.6. For 
blade-passage frequency at the 3 speeds 9,600, 1 1,400, and 12,000 a least squares fit of the micro- 
phone data to 2 modes was used. The harmonic, 2BPF, at 9,600, having 3 propagating modes could 
be resolved into these 3 modes. 


Table 2 gives the inputs, p M , and the fan source mode structure, Pp, for the selected rig speeds. At 
9600 rpm, the BPF mode ju = 1 is cut off. For this reason there is a larger difference between p M and 
Pp for this case than for the remaining modes which are all cut on. The fan face level, pp , is also 
higher than the level at the microphone, p M . Well cuton modes show smaller differences between 
Pp and p M with the p M levels larger. This results from the variation of geometry and flow velocity 
throughout the inlet. Figures 13 and 14 show examples of mode amplitude variations within the 
duct. Modes just above cutoff, such as BPF mode /i=l at 11,400 rpm and 2BPF mode fx=2 at 9,600 
rpm, result in larger differences between pp and p M than well cut on modes. 


If, in future programs, it is deemed essential to obtain information about decaying modes generated 
by the fan, further sensitivity studies must be made. For example, it may turn out that a different 
inlet configuration, providing a closer fan-to-microphone spacing (lower value of x/r) will work. 
The so-called “short inlet” should give a value of x/r of about 0.25. Computer simulation studies 
could be made to determine if reliable results can be extended to 4 modes. 
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CONCLUSIONS AND RECOMMENDATIONS 


In using measurements from rotating microphones in fan inlets to predict far field directivity pat- 
terns, waves reflected from the inlet mouth back to the microphones can lead to significant errors 
in the predictions. A procedure has been derived for minimizing these effects by use of a finite ele- 
ment radiation code that separates upstream-going waves from downstream-going waves. 

Via an analytical simulation, it was shown, for a mode near cutoff, that far field projections could 
vary erroneously by large amounts, depending on where the rotating microphone plane happened 
to be in the internal standing wave pattern. The enhanced data reduction scheme successfully col- 
lapsed the error band to about 0.5 dB. The system was then applied to four test runs of the 17 inch 
ADP fan model where the rotating mic data were available. Comparison of far field projections us- 
ing the basic and enhanced schemes yielded nearly the same results, apparently because the mic array 
was not near a node in the standing wave pattern. For conditions that include modes near cutoff and 
mic array locations near standing wave nodes, the enhanced scheme can make a significant improve- 
ment in accuracy of far field projections. A criterion is needed for when the enhanced scheme is 
needed (based, say on cutoff ratio). 


It is also recommended that a modification to the radiation code be examined for feasibility: The 
present code uses a common plane as the source of both the noise and flow and the same grid for 
each field. Studies such as those discussed in this report require that the location of the noise source 
be taken at various axial locations within the duct The result is that the details of the flow field 
change slightly as the location of the noise input is changed. To eliminate this, means should be 
found to allow separate input planes for the noise and flow, possibly by interpolating the flow field 
onto the desired noise grid. 

For the Application 2 task, which consists of determining the fan mode structure, pp , from micro- 
phone measurements, p M , using measured data where the normalized fan to microphone distance, 
x/r, was 0.6, it is concluded that only propagating or weakly decaying modes be employed. 

If in other configurations, the fan and microphone planes are closer together (x/r <0.2) it is con- 
cluded that decaying modes can also be found. 

It is recommended that, in planning other tests and fan configurations, sensitivity studies be made 
at an early stage to determine what features are needed to allow recovery of the desired information. 

A correlary recommendation, applicable to existing test data on the 17” ADP rig, is to analyze the 
so-called short inlet configuration to see if the closer fan-to-microphone spacing will allow recov- 
ery of the decaying modes. 
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APPENDIX A 
NOTATION 


a,b,c,d 
ai ••• &6 
BPF 
Cmu 
D,D- ! 

Dmf’ 

expi(z) 

I 

kmn 

m 

P 

P* 

P 

Pf 

Pm 

Pm 

Pm 

P 

Pi. P2» P3 
Pl_’P2_’V? 

Pi . P2 , P3 

Pp 

Qm|i 

r 

t 

T 

x 


constants 

coefficients of interpolation polynomial 
fan blade passage frequency (Hz.) 

coefficient of (m,/c) circumferential-radial mode 
calibration matrix, inverse 

calibration matrix relating p at mic plane to p + at fan face 
calibration matrix relating p at mic plane to p + at mic plane 
complex exponential function, e a 
unit matrix 

Bessel function of the first kind of order m 

eigenvalue for (m,/t) mode 

circumferential mode index 

total acoustic pressure, p + + p “ 

upstream traveling acoustic pressure component 

downstream traveling acoustic pressure component 

upstream traveling acoustic pressure component at fan face 

total acoustic pressure at mic plane 

upstream traveling acoustic pressure component at mic plane 
downstream traveling acoustic pressure component at mic plane 
complex acoustic pressure 

mode components sum of upstream and downstream components 

upstream traveling mode components 

downstream traveling mode components 

mode amplitude coefficient = Cmu 

eigenvalue for (m ,/u) mode 

radial coordinate 

time 

vector transpose (superscript) 
axial coordinate 

Bessel function of the second kind of order m 


0 

t* 


'lV(r) 


co 


circumferential coordinate 
radial mode index (0,1, ...) 
mode radial shape function 
radian frequency = 2jc x frequency 
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APPENDIX B 
SENSITIVITY STUDIES 


The general basis for these studies is to use computer simulation with postulated, known fan modes 
pp as inputs to the radiation code in order to obtain the outputs, p M . These outputs were perturbed 
by small amounts and the perturbed quantities were inserted into Equation 13 to simulate the prob- 
lem of obtaining pp from p M . The resulting pp values and their deviations from the known input 
provided a guide to the effects of small variations in the “data”. Rather than using just the test micro- 
phone location at x/r = 0.6, the process was applied over the range of mic. locations between x/r = 
0 to x/r = 0.8 in steps of 0.05 to find how the results varied with source-to-microphone separation. 

At each rig speed one of a sequence of unit pp modes, [1,0,0,0] T , [0,1,0,0] T , etc. was applied. 
Values of p M were obtained from the radiation code, and were perturbed according to the following 
scheme: 

It was assumed that microphone reading errors could be expressed as a fraction of 1 decibel and these 
were applied to the modal coefficients, p M , as follows. 


If lp M l = 

a 

(Pascals) and the perturbations were arbitrarily taken as 

0.5 


b 


-0.4 


c 


0.3 


_d_ 


0.1 


then the values of p M would be replaced by multiplying them by factors of 10 dB / 20 , giving the per- 
turbed vector 


10-V 20 a 

— 

1.059 a 

10 ~ - 4 / 20 b 


0.955 b 

10 3 /20 c 


1.035 c 

-10 1 / 20 d- 


- 1.012 d- 


(Pascals) 


(Note a,b,c,d are the amplitudes of the complex quantities; the phase angles were left unchanged.) 

Table 3A gives the results of the process for BPF at 12,000 rpm. The fan mode exciting the system 
is a unit magnitude pi = 0 mode, [1,0,0,0] T . If the resulting values of p M were not perturbed, the 
computed values of pp would be exactly equal to the input, [1,0 ,0,0] Pascals. In decibel units, the 
tabulated values in the successive colums would be 0 dB, - oo dB, - « dB, and - oo dB, and these 
would hold for all distances x/r from 0 to 0.8. 

With the perturbed values of p M , examination of the first 2 or 3 lines, where the microphone and 
fan face are close together reveals two features: The first column values, corresponding to the pi - 0 
unit input mode are very close to the assumed perturbation of 0.5 dB. The last 3 column values are 
at least 40 dB down, which is very good from a practical viewpoint. 
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Switching attention now to Table 3B, which gives results for unit // = 1 mode, [0,1,0, 0] T , the second 
column values are close to the perturbation of -0.4 dB : all other columns are way down for x/r small. 
Corresponding results are obtained for unit n = 2 and 3 inputs, as shown in Tables 3C and 3D. 


In summary, for sufficiently small x/r, the process of obtaining the fan mode structure pp from mi- 
crophone readings, p M subject to small error, is very well-behaved; errors in the result are propor- 
tionate to errors in the measured values p M . 

However, as x/r increases, the nature of the results changes markedly - calculated values of all modes 
grow, until at the upper limit, x/r = 0.8, the listed pressures of all modes are manifestly absurd. This 
progressive degeneration holds for all 4 excitations, shown in Figures 3A thru 3D. For the micro- 
phone location x/r = 0.6 where actual rig data were acquired, examination of the set of Figures 2 
shows obviously unacceptable results. 

Similar results were obtained for BPF at 9,600 rpm and 11,400 rpm, and for 2BPF at 9,600 rpm. 
In addition, the entire set of calculations was repeated using different perturbations: 0.3dB, 0. ldB, 
0.5dB and -0.4dB. These showed the same trends as just described; good results for small x/r, de- 
grading to very poor information at the test microphone location, x/r = 0.6. 

It is clear that any future testing of fan configurations should be preceded by computer-simulated 
runs using the appropriate parameters in the radiation code. By this means, the test microphone loca- 
tion can be selected in advance to give reliable predictions of the fan source strength, pp . 

In an attempt to find some useful information about source excitation, pp , from the ADP rig data 
on hand, it was conjectured that if resolution of the microphone data by least squares fit to 4 modes 
were to be fit instead to fewer, lower order modes, the results might be a valid measure of the propa- 
gating modes present. It may be that decaying modes at x/r = 0.6 are sufficiently contaminated by 
system “noise” so as to give a false prediction of the source, and that by simply eliminating them 
from the model, the remaining modes could be found reliably. Fortunately, this turned out to be true, 
as shown in Tables 4,5,6, and 7. 

Table 4 applies to BPF at 9,600 rpm. Only the n = 0 mode propagates, but the n - 1 mode does 
not decay very rapidly and was included in a least squares fit. The perturbations were 0.5 dB and 
-0.4 dB, as before. (Modes /r = 2 and n = 3 and their perturbations are simply absent throughout 
the process. Table 4A gives results vs. x/r for unit fi = 0 input, and Table 4B is for unit p - 1 mode. 
Results are clearly satisfactory at x/r = 0.6 (The results initially obtained for a 4-mode fit at this speed 
were useless at x/r = 0.6; for brevity they are not shown in this report.) 

Tables 5 A and 5B are for 11,400 rpm, where the first two modes propagate. 

Tables 6Aand 6B represent results at 12,000 rpm, where there are still 2 propagating modes. Again 
results at x/r = 0.6 are satisfactory. 
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Tables 7 A, B, C apply to twice blade passage frequency at 9,600 rpm. Here, the first 3 modes are 
propagating and a corresponding 3-mode least squares fit was used for p M . Three tables, A, B, C, 
are now needed to display the result of perturbations 0.5, -0.4, and 0.3. It is evident that results are 
good, not only at x/r = 0, but all the way to 0.8. 

These sensitivity studies also reveal information about Application 1 type tests. These tests involve 
extracting the forward component, pjJ at the microphone plane, rather than at the fan, pp , from 
total microphone pressure, p M . The sensitivity studies show that when the plane, p M is close to the 
p + plane the prediction scheme using all 4 modes works well. In Application 1, these planes are 
coincident (x/r ~ 0) so that the results obtained there can be regarded with confidence. 
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TABLE 1 


Fan— generated Mode Amplitudes, pp , Computed from Measured Microphone Values, p^, 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Blade Passage Frequency 


9,600 rpm 11,400 rpm 12,000 rpm 


Mode No. 

Pm 

Pf 

Pm 

Pf 

Pm 

Pf 

0 

33.5 

62.2 

31.6 

62.1 

28.4 

69.1 

1 

24.4 

63.3 

27.8 

66.5 

27.0 

73.2 

2 

16.6 

79.3 

13.6 

76.7 

35.8 

84.5 

3 

17.6 

99.4 

20.8 

98.2 

31.2 

104.8 
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TABLE 2 


Fan Modes, pp, from variable Least Squares fit of Microphone Modes, p M 
All Amplitudes have Units of Decibels re. One Pascal. 


BPF 

9,600 rpm 


BPF 

11,400 rpm 


BPF 

12,000 rpm 


2BPF 
9,600 rpm 


Mode No. 


/*• 

Pm 

Pf 

Pm 

0 

33.5 

30.4 

31.7 

1 

24.0 

41.2 

27.2 


2 


Pf 

Pm 

Pf 

Pm 

Pf 

28.1 

27.9 

24.1 

47.1 

47.2 

22.5 

25.9 

24.9 

45.0 

42.7 




48.3 

40.3 
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TABLE 3A 


Fan-generated Mode Amplitudes, Pp , Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Blade Passage Frequency 
12,000 rpm 

Input mode, p^ = [1,0,0,0] T 


Distance x/r 


0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 


Mode /<= 0 


0.4983 

0.4968 

0.4993 

0.4968 

0.4862 

0.4620 

0.4428 

0.4755 

0.6716 

1.1256 

1.4583 

2.7357 

0.3406 

11.8982 

21.2753 

28.7842 

39.5239 


Mode 1 


-42.0939 

-42.1069 

-44.0812 

-52.6600 

-39.9090 

-33.5101 

-32.0718 

-33.1389 

-27.6733 

-18.7334 

-12.7149 

-3.0345 

-0.2393 

14.2791 

24.4680 

32.5388 

43.5666 


Mode /a= 2 


-60.6872 

-51.3050 

-45.1004 

-42.4603 

-53.1564 

-35.5056 

-27.3077 

-22.6374 

-15.4206 

-7.2096 

-0.6964 

8.5145 

11.4179 

26.1699 

36.3513 

44.4155 

55.4345 


Mode /t=3 


-63.9803 

-51.5325 

-42.5006 

-34.7413 

-26.5955 

-19.5173 

- 11.2002 

-2.5900 

6.1875 

14.2228 

20.3841 

28.9615 

31.2333 

46.0359 

56.2520 

64.3228 

75.3530 


0=0, n=\ propagating) 



TABLE 3B 


Fan-generated Mode Amplitudes, pp. Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Blade Passage Frequency 
12,000 rpm 

Input mode, pp = [0,1,0,0] T 


Distance x/r 


0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 


Mode fi=0 


-60.7032 

-47.6168 

-42.7402 

-42.6137 

-42.6407 

-40.5633 

-37.3929 

-29.7582 

-24.1306 

-18.2286 

-11.5820 

6.7835 

6.9592 

15.6758 

23.8776 

32.2267 

43.1881 


Mode fi = 1 


-0.4001 

-0.3984 

-0.3986 

-0.3851 

-0.3419 

-0.2695 

-0.2245 

-0.3975 

-0.9456 

-1.7903 

-0.3011 

10.0674 

12.8882 

20.5211 

28.2255 

36.3542 

47.2124 


Mode ju=2 


-42.3314 

-36.9705 

-32.0085 

-27.1599 

-21.8361 

-17.4699 

-13.7425 

-10.1966 

-6.5995 

-1.3621 

5.5007 

22.0657 

23.2836 

31.8811 

39.9511 

48.1713 

59.0526 


Mode n= 3 


-56.9170 

-46.7569 

-37.5835 

-29.2593 

-21.0339 

-13.4298 

-5.6653 

1.6835 

9.1528 

16.9708 

24.4950 

42.6680 

42.8051 

51.4782 

59.6682 

68.0105 

78.9717 


0=0, ju=l propagating) 



TABLE 3C 


Fan-generated Mode Amplitudes, pp , Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Blade Passage Frequency 
12,000 rpm 

Input mode, pp = [0,0,1,0] T 


Distance x/r 


0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 


Mode fi = 0 


-71.9642 

-70.3679 

-68.8556 

-68.6473 

-65.3853 

-61.1792 

-61.8010 

-60.2849 

-52.8500 

-42.9073 

-34.2576 

-17.4971 

-16.2695 

-7.3879 

1.3164 

9.6988 

20.2520 


Mode ju=l 


-60.2117 

-57.6402 

-54.2661 

-49.9164 

-46.1812 

-45.3530 

-46.3129 

-46.1968 

-44.4384 

-38.5787 

-30.5533 

-13.6060 

-12.3818 

-3.4279 

5.3137 

13.7048 

24.2569 


Mode n=2 


0.2996 

0.2998 

0.2986 

0.2933 

0.2767 

0.2431 

0.2044 

0.2028 

0.3137 

0.5832 

0.8524 

4.8727 

2.3403 

7.6714 

16.4532 

25.1665 

35.9951 


Mode [i = 3 


-64.7547 

-59.2213 

-52.7881 

-46.2554 

-40.4471 

-36.6649 

-33.8795 

-26.1994 

-15.8930 

-6.3098 

1.9737 

18.3904 

19.5278 

28.3869 

37.0983 

45.4828 

56.0359 


(/t=0, /i = 1 propagating) 



TABLE 3D 


Fan-generated Mode Amplitudes, pp , Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Blade Passage Frequency 
12,000 rpm 

Input mode, pp = [0,0,0, 1] T 


Distance x/r 


0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 


Mode /i = 0 


-66.6672 

-74.0001 

-68.0009 

-65.1911 

-64.5384 

-65.2605 

-67.5697 

-69.9456 

-71.8641 

-59.9777 

-49.2954 

-25.8315 

-26.0609 

-15.6523 

-6.1688 

2.2951 

12.5441 


Mode ju=l 


-66.3373 

-61.4967 

-59.8443 

-59.4957 

-59.1168 

-59.2930 

-59.5946 

-56.9243 

-53.6292 

-50.3142 

-44.9420 

-21.8919 

-22.1812 

-11.6949 

-2.1730 

6.3003 

16.5493 


Mode fi= 2 


-57.5148 

-56.1559 

-54.0236 

-51.6988 

-48.9465 

-46.0591 

-43.8164 

-42.6455 

-42.6046 

-40.4912 

-31.7806 

-10.4258 

-9.6313 

0.5281 

9.8458 

18.2210 

28.4219 


Mode n = 3 


0.1036 

0.1032 

0.1021 

0.0973 

0.0846 

0.0749 

0.0920 

0.2276 

0.4047 

0.3966 

-0.6079 

10.6707 

6.7499 

19.1797 

29.3193 

37.9985 

48.3259 


(ju=0, ju=l propagating) 



TABLE 4A 


Fan-generated Mode Amplitudes, pp , Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Blade Passage Frequency 
9,600 rpm 

Input mode, Pp = [1,0] T 


Distance x/r 


0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 


(ju=0 propagating) 


Mode ju=0 

Mode /u=l 

0.4979 

-39.5877 

0.4975 

-36.7478 

0.4962 

-34.1682 

0.4935 

-31.7672 

0.4874 

-29.3507 

0.4780 

-27.2583 

0.4682 

-25.5668 

0.4630 

-24.7413 

0.4649 

-24.3966 

0.4772 

-24.5669 

0.5000 

-24.5604 

0.5484 

-21.3420 

0.5361 

-20.2751 

0.5051 

-18.2635 

0.4331 

-15.2950 

0.3703 

- 7.6076 

0.5021 

1.3457 
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TABLE 4B 


Fan-generated Mode Amplitudes, Pp , Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Blade Passage Frequency 
9,600 rpm 

Input mode, pp = [0,1] T 


Distance x/r 


0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 


Mode ji=0 


-51.9017 

-50.3314 

-49.6326 

-49.1280 

-47.0167 

-44.3997 

-43.0632 

-42.7220 

-43.1054 

-43.6889 

-44.2512 

-40.8661 

-44.7555 

-45.3858 

-46.2091 

-46.8506 

-44.4952 


Mode fi = 1 


-0.3989 

-0.3984 

-0.3970 

-0.3940 

-0.3873 

-0.3768 

-0.3661 

-0.3604 

-0.3625 

-0.3760 

-0.4011 

-0.4552 

-0.4414 

-0.4066 

-0.3275 

-0.2578 

-0.3723 


(iu = 0 propagating) 



TABLE 5A 


Fan-generated Mode Amplitudes, pp , Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Blade Passage Frequency 
11,400 rpm 

Input mode, pp = [1,0] T 


Distance x/r 


0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 


Mode n=0 


0.4958 

0.4951 

0.4937 

0.4888 

0.4827 

0.4802 

0.4789 

0.4647 

0.4391 

0.4336 

0.4506 

0.4846 

0.4931 

0.4943 

0.4860 

0.4708 

0.4384 


Mode fx-\ 


-36.3481 

-38.1990 

-38.0061 

-37.0658 

-35.8805 

-35.3279 

-35.9649 

-31.9297 

-28.3415 

-28.5792 

-30.3333 

-34.0732 

—44.9382 

-40.1698 

-33.6664 

-31.4760 

-30.9866 


0=0, n - 1 propagating) 



TABLE 5B 


Fan-generated Mode Amplitudes, pp\ Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Blade Passage Frequency 
11,400 rpm 

Input mode, pp = [0,1] T 


Distance x/r 


0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 


Mode n = 0 


-51.9393 

-50.8032 

-47.4501 

-41.7832 

-37.7854 

-35.9670 

-35.8325 

-35.5587 

-34.6581 

-33.8977 

-34.6274 

-28.7592 

-38.3118 

-36.1119 

-33.1425 

-31.7557 

-30.8946 


Mode 1 


-0.3965 

-0.3959 

-0.3943 

-0.3889 

-0.3821 

-0.3794 

-0.3778 

-0.3623 

-0.3342 

-0.3282 

-0.3468 

-0.3837 

-0.3936 

-0.3949 

-0.3855 

-0.3686 

-0.3333 


(/z=0, fi=\ propagating) 



TABLE 6A 


Fan-generated Mode Amplitudes, Pp, Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Values have Units of Decibels re. One Pascal. 


Blade Passage Frequency 
12,000 rpm 

Input mode, pp = [1,0] T 


Distance x/r 

Mode fi=Q 

Mode //=1 

0.00 

0.4983 

-41.5396 

0.05 

0.4969 

-41.3311 

0.10 

0.5001 

-43.2448 

0.15 

0.4991 

-54.1553 

0.20 

0.4914 

-40.6982 

0.25 

0.4733 

-33.7336 

0.30 

0.4560 

-31.3059 

0.35 

0.4559 

-31.6531 

0.40 

0.4618 

-33.0684 

0.45 

0.4628 

-32.1487 

0.50 

0.4552 

-28.9901 

0.55 

0.5056 

-28.1050 

0.60 

0.4890 

-34.2181 

0.65 

0.4989 

-47.0594 

0.70 

0.4806 

-33.4832 

0.75 

0.4516 

-29.8826 

0.80 

0.4061 

-29.0436 


0=0, n=\ propagating) 
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TABLE 6B 


Fan-generated Mode Amplitudes, pp, Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Values have Units of Decibels re. One Pascal. 


Blade Passage Frequency 
12,000 rpm 

Input mode, pp = [0,1] T 


Distance x/r 


0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 


Mode fi=0 


-60.3535 

-47.7771 

-42.4365 

-41.4934 

-40.3622 

-36.6253 

-34.2635 

-33.3539 

-33.8858 

-34.3899 

-34.6144 

-32.5263 

-41.6366 

-43.9171 

-34.7168 

-31.0264 

-29.1392 


Mode /u=l 


-0.3993 

-0.3978 

-0.4014 

-0.4003 

-0.3917 

-0.3717 

-0.3526 

-0.3525 

-0.3590 

-0.3601 

-0.3516 

-0.4067 

-0.3891 

-0.4000 

-0.3797 

-0.3474 

-0.2978 


0=0, n=\ propagating) 



TABLE 7A 


Fan-generated Mode Amplitudes, Pp , Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Twice Blade Passage Frequency 
9,600 rpm 

Input mode, pp = [1,0,0] T 


Distance x/r 

Mode fi=Q 

Mode ju=l 

Mode ju=2 

0.00 

0.4981 

- 57.6101 

-60.5270 

0.05 

0.4972 

-54.2157 

-60.7388 

0.10 

0.4963 

-50.0314 

-69.3182 

0.15 

0.4931 

-43.0440 

-56.8064 

0.20 

0.4842 

-39.2625 

-50.6235 

0.25 

0.4782 

-37.9375 

-55.2916 

0.30 

0.4701 

-35.7663 

-51.9492 

0.35 

0.4541 

-33.6034 

^42.3256 

0.40 

0.4453 

-33.0857 

^40.4393 

0.45 

0.4419 

-32.4335 

-42.4967 

0.50 

0.4395 

-31.5047 

-41.2395 

0.55 

0.4504 

-32.2531 

-40.7099 

0.60 

0.4755 

-34.7792 

-43.6341 

0.65 

0.4932 

-38.6193 

-49.2649 

0.70 

0.4941 

-66.9818 

-48.2364 

0.75 

0.4566 

-35.4951 

-38.1099 

0.80 

0.3720 

-30.7390 

-33.2868 


(ju= 0, fi= 1, ju=2 propagating) 
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TABLE 7B 


Fan-generated Mode Amplitudes, pp , Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Twice Blade Passage Frequency 
9,600 rpm 

Input mode, pp" = [0,1, 0] T 


Distance x/r 


0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 


Mode n=0 


-59.7327 

-48.4808 

-44.2941 

-42.0531 

-37.1942 

-35.2459 

-34.4992 

-33.2442 

-32.7316 

-32.6315 

-32.5768 

-33.3115 

-36.0627 

-39.3813 

-39.6106 

-33.5685 

-29.8096 


Mode n=\ 


-0.3981 

-0.3966 

-0.3922 

-0.3777 

-0.3689 

-0.3510 

-0.3116 

-0.2498 

- 0.2221 

-0.2310 

-0.2427 

-0.2681 

-0.3199 

-0.3613 

-0.3492 

-0.2667 

-0.1215 


Mode [i=2 


-44.6300 

-42.4633 

-38.2509 

-36.7339 

-37.7396 

-35.1722 

-31.2076 

-29.7316 

-28.9187 

-29.0618 

-30.0752 

-31.2239 

-32.5301 

-33.5881 

-32.2830 

-30.7972 

-30.2460 


(/t=0, [i=l, [i=2 propagating) 



TABLE 7C 


Fan-generated Mode Amplitudes, pp , Computed from Perturbed Microphone Values, p M , 
Using Equation 13. All Amplitudes have Units of Decibels re. One Pascal. 


Twice Blade Passage Frequency 
9,600 rpm 

Input mode, Pp' = [0,0, 1] T 


Distance x/r 


0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 


Mode /i=0 


-72.6975 

-63.0772 

-60.5412 

-64.4835 

-64.4049 

-58.2798 

-50.5067 

-46.2464 

-45.6158 

-45.2430 

-44.8579 

-46.6272 

-53.3277 

-67.2021 

-47.1183 

-39.1491 

-34.2684 


Mode fi=l 


-50.7547 

-53.0153 

-48.2286 

-39.8584 

-36.9821 

-36.7788 

-35.3777 

-32.8241 

-32.3553 

-32.9497 

-33.3355 

-33.3128 

-34.7917 

-37.4127 

-35.4580 

-32.6590 

-31.8111 


Mode fi = 2 


0.2979 

0.2975 

0.2942 

0.2840 

0.2850 

0.2746 

0.2461 

0.2042 

0.1869 

0.1990 

0.2125 

0.2253 

0.2483 

0.2687 

0.2566 

0.2176 

0.1657 


(H=Q, ju=l, ju=2 propagating) 
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Figure 1. Comparison of far field directivities for p and p + inputs at four axial locations. 

BPF at 7500 rpm. Mode = 6,0. Noise field obtained from unit input amplitude at 
the fan face. 
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Measure 



Pm - Pm + Pm 


Backward Wave Component 


Required Forward Wave Component 


Relation p M = D M mPm + determined by “calibration” with inlet radiation 
code. 

In expanded form (for 3 radial modes): 

Complex Pi] = Du D 12 D 13 ] # T pi + " 

radial mode p 2 D 21 D 22 D 23 p 2 + 

vector p3_ m _ D 31 D 32 D 33 _ mm _ p 3 + _ m 

When elements of D M m are known, any measured p M can be resolved, 

Pm + = Dmm -1 Pm 


Figure 2. Application 1 - Method of obtaining forward waves, p^, at the mic array from 
measured total pressures, p^, at the same plane. 
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Figure 3. 
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Mesh system used to obtain forward waves, p^, at the mic array from measured total 
pressures, p^, at the same plane. 
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Pm - D mm Pm + 
write as p = D p + 


for 3 x 3 

Pi 

= 

' Du 

D12 

D13 


Pl + 


P 2 


D21 

D22 

D23 


P 2 + 


P 3 


D31 

D32 

D33 


P 3 + 


Calibrate with successive known, unit modal inputs 


first 

D11" 

— 


D21 



P 31 _ 



second 

D12 

= 


D 2 2 



D32 

- 

third 

D13 

_ 


D23 



D33 




First unit input gives 
first column elements 
of D 


Second unit input gives 
second column of D 


Third unit input gives 
third column of D 


Outputs from computer program are columns of Dmm matrix 


Figure 4. Computer calibrations to give elements of D^. 
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FACE MIC 

FACE PLANE 


(required) Pp 

Pf" 



Pm 



Relation p M = D MF p+ 


In expanded form (for 3 radial modes): 


Complex 

’ Pi" 

= 

Du D -12 Dl 3 

• 

Pl + 

radial mode 

P2 


D21 D22 D23 


P2 + 

vector 

- P 3 - 

M 

CO 

CO 

Q 

C\J 

CO 

Q 

T“ 

CO 

Q 
1 

MF 

L p 3 + J 


When elements of D MM are known, any measured p M can be resolved, 

Pm + = D mf _1 Pm 


Figure 5. Application 2 - Method of obtaining forward waves, pp , at the fan face from 
measured total pressures, p M , at the plane of the rotating microphones. 
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Figure 6. Mesh system used to obtain forward waves, pp , at the fan face from measured 
total pressures, p M , at the plane of the rotating microphones. 
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Pm = Dmf Pf + 
write as p = D p + 


for 3x3 

Pi 

= 

" Du 

D12 

D13 


P 1 + 


P 2 


D 2 i 

D22 

D23 


P 2 + 


P 3 


D 3 i 

D32 

D33 


P 3 + 


Calibrate with successive known, unit modal inputs 


first 

D11" 

=5 


D 2 i 



_ D 31 _ 



second 

D12 

= 


D22 



D32 

- 

third 

D 13 




D 23 



_D 33 _ 




First unit input gives 
first column elements 
of D 


Second unit input gives 
second column of D 


Third unit input gives 
third column of D 


Outputs from computer program are columns of Dmf matrix. 


Figure 7. Computer calibrations to give elements of D^. 
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Figure 10. Estimated far-field directivity from measured rotating microphone data. 
BPF, 11,400 rpm corrected rotor speed, circumferential mode = -6. 
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Figure 11. Estimated far-field directivity from measured rotating microphone data. 
BPF, 12,000 rpm corrected rotor speed, circumferential mode = -6. 
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Figure 12. Estimated far-field directivity from measured rotating microphone data. 

2BPF, 9,600 rpm corrected rotor speed, circumferential mode =10. 
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NORMALIZED AXIAL LOCATION, X 

Figure 13a. Comparison of mode amplitudes, p, at various axial locations 
for unit input p + modes. BPF = 9,600 rpm. Input p + = 0. 
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P = 2 
P = 3 





NORMALIZED AXIAL LOCATION, X 

Figure 1 13b. Comparison of mode amplitudes, p, at various axial locations 
for unit input p + modes. BPF = 9,600 rpm. Input p, + = 1. 


r 




SOUND PRESSURE LEVEL -dB 


NORMALIZED AXIAL LOCATION, X 

Figure 13c. Comparison of mode amplitudes, p, at various axial locations 
for unit input p + modes. BPF = 9,600 rpm. Input p + = 2. 


P= 2 <^, 



Kg a 




0.0 0.1 


O.S 0.6 


NORMALIZED AXIAL LOCATION, X 


Figure" 13d. Comparison of mode amplitudes, p, at various axial locations 
for unit input p + modes. BPF = 9,600 rpm. Input p 3 + = 1. 





Figure 14a. Comparison of mode amplitudes, p, at various axial locations 



NORMALIZED AXIAL LOCATION, X 


Figure 14b. Comparison of mode amplitudes, p, at various axial locations 
for unit input p + modes. BPF = 12,000 rpm. Input p* = 1. 
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NORMALIZED AXIAL LOCATION, X 

Figure 14c. Comparison of mode amplitudes, p, at various axial locations 



Figures 14d. Comparison of mode amplitudes, p, at various axial locations 
for unit input p + modes. BPF = 12,000 rpm. Input p 3 + = 1. 
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